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Abstract 

A simple mathematical expression for the universal map for cellular automata is 
found in closed form with the help of a digit function, whose most basic properties 
are established. This result is found after proving a theorem on the composition of 
functions on finite sets. The expression (and the technique used to obtain it) opens 
the possibility of gaining mathematical insight in any cellular automaton rule since 
it constitutes at the same time a simple and fast algorithm to implement any such 
rule. 
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1 Introduction 

Cellular automata (CA) constitute an important tool in the study of complex systems [T1I21I31H] . 

CA have been satisfactorily used to model physical [5], chemical [6] and biological processes and, 
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in certain cases, constitute an alternative to continuous models described by partial differential 
equations [7j. Natural patterns as those found in some seashells, like the Conus marmorea |H] can 
readily be obtained with CA [1] . Although such patterns may involve many microscopic degrees 
of freedom, a coarse graining of them can also be generally described by cellular automata [9llT0] 
and the latter can thus be used to capture the collective behavior of spatiotemporal pattern 
forming systems. 

Despite involving only a hnite number p of dynamical states and the interactions having a hnite 
range, already the most elementary CA can generate the highest complexity mM- The mathe¬ 
matical description of such systems in terms of discrete maps seems to have been only recently 
attempted mn] and previous researchers have mostly relied upon the computer to get insight 
in CA behavior. In the monumental work [1] , Chua has pointed out how useful is to have math¬ 
ematical expressions to gain insight on CA behavior: nonlinear dynamical methods, number 
theoretical, algebraic and statistical techniques can then be fruitfully applied to these systems. 
Chua describes a universal neuron which provides a mathematical description for Wolfram’s 
256 elementary CA |1]. However, this model generally depends on eight adjustable parameters 
that need to be sought with the computer, and which are then tabulated occupying several 
pages in [1]. Ozelim et. ah m have also introduced the so-called iota-delta function to describe 
Wolfram’s elementary CA. Again, however, their model depends on freely adjustable param¬ 
eters that need to be tabulated separately: in order to describe CA behavior, such numerical 
parameter values need to be given as input together with the specihcation of the dynamical 
state of the cells. In our opinion, the presence of such parameters in the equations is an obvious 
major drawback from a mathematical point of view. 

In [I2] we have introduced a universal map for cellular automata which does not depend on any 
adjustable parameter and which is not only able to simulate Wolfram’s 256 elementary CA, 
but any deterministic CA with any number of dynamical states. This map can also be easily 
extended to any number of dimensions. In |I3] and HI we have discussed some symmetries 
of the universal CA map that help to classify the rules into symmetry classes |T3] and to 
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understand the kind of symmetry breaking that leads to complexity at this elementary level 
[T^ . We have termed i3-calculus the mathematical techniques introduced in [T2] because these 
strongly rely on a function of two arguments that is called the boxcar or i3-function [12] 


B{x,y) 


1 ( X + y 

2 V|a; + 2/| 


x-y 

x-y\ 


( 1 ) 


This function constitutes the building block of i3-calculus. It returns the sign of the second 
argument y, if —12/| < x < \y\, sign y/2 if x = y and ’0’ otherwise. The S-function is thus 
an even function of its hrst argument and an odd function of its second one. The function is 
represented in Fig. [Tj 


As introduced in na, the universal map is given in terms of a sum that involves terms, 

where p is the alphabet size and I and r denote the number of neighbors to the left and to the 
right (respectively) of the symbol that is updated at the next time step. In this paper we reduce 
this sum to a single evaluation over the neighborhood values: the universal map is fully reduced 
to a closed form with mathematical structure. The approach makes use of a digit function whose 
main mathematical properties are hrst established. Then, we show how the digit function alone 
is able to describe all deterministic CA in ID. Extensions to more dimensions also become 
trivial although we refrain from discussing them here. 


We briehy describe the universal map for cellular automata, whose closed form we then derive. 
Let us consider a ID ring containing a total number of Ng sites. An input is given as initial 
condition in the form of a vector xq = (xj,..., x^®). Each of the Xg is an integer in [0,p — 1] where 
superindex j G [1, W] specihes the position of the site on the ID ring. At each t the vector 
Xi = (xj, ...,x^") specihes the state of the CA. Periodic boundary conditions are considered so 
that xf^'*"'"^ = x] and x° = x^‘. Let be taken to denote the value of site j at time step 
t + 1. Formally, its dependence on the values at the previous time step is given through the 
mapping xi_^_i = ’‘Rp{x{'^\ ■■■Xt, ■■■,xi~^), which we abbreviate as xi_^_i = ’‘Rp{x{) with the 
understanding that the function on the r.h.s depends on all site values within the neighborhood, 

with range p = / + r + 1, which contains the site j updated at the next time {I and r denote the 
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number of cells to the left and to the right of site j respectively). We take the convention that 
j increases to the left. The integer number n in base 10, which runs between 0 and 
indexes all possible neighborhood values coming from the different conhgurations of site values. 
Each of these conhgurations compares to the dynamical conhguration reached by site j and its 
r and I hrst-neighbors at time t and given by 


ni = ( 2 ) 

k=—r 


We will refer to this latter quantity often as the neighborhood value. The possible outputs a„ 
for each conhguration n are also integers G [0,p — 1]. An integer number R can then be given 
in the decimal base to fully specify the rule ’‘R^ as 


pT + l + l_l 

R= ( 3 ) 

n=0 


This is the so-called Wolfram code of the CA dynamics. With all these specihcations we have 
the following universal map [12] 


„r+Z + l 


-1 


xi+i = ^Rp{xl) = n - 


n=0 


I 

E 

k=—r 


„fc+r i+fc 

P Xt , 


= CL^ 


( 4 ) 


This map, although describing all ID deterministic cellular automata which are hrst order in 
time has the shortcoming that it generally depends on a sum which can contain up to 
terms. Although Eq. Q is advantageous in a certain sense because the ^-functions in the map 
behave as an orthonormal base (thus making the transition to a Hilbert space clear, if the a„ 
are allowed to take on complex values) the map is computationally demanding for / r -|- 1 of 
p large and a closed form can be desirable in certain cases to gain further mathematical insight 
and to test the validity of certain calculations with the computer. The Ending of such closed 
form is the goal of the next section. 
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2 The digit function 


Let now p > 1, p eN. The representation in radix p of a; G M is a (generally infinite) sequence 
of the form 

X = unP^ + + ... + ao + a-ip“^ + a-2P~^ + • • • (5) 

where the Ofeg^’s, the so-called digits of the expansion, are all positive or negative, with absolute 
values being integers in S, the set of integer numbers between 0 and p — 1. Every real number 
can be represented in this form. Let [.. .J denote the closest lower integer (fioor function) and 
[...] the closest upper integer (ceiling function). The following function 

sign(a:)p“'^[p^|a:|J = a^p^ + aAr_ip^“^ + ■■ ■ + ao + a-ip“^ + a_ 2 P~^ a_DP~^ (6) 


is the rational truncation to precision D in radix p of the real number x. We thus realize the 
following interesting fact: the operators [.. .J and multipliying by p^ do not commute |T5], i.e. 
.. .J 7 ^ p-^ [.. .J. Indeed, we have 

sign(x) {p~^[p^\x\\ - L|x|j) = a_ip“^ a_ 2 p“^ a-DP~^ (7) 


which gives the rational truncation to precision D in radix p of the fractional part of x. Let us 
consider x nonnegative for simplicity. Now, note that 

= aNP^~^ -t- Ofc+ip Ofc (8) 

and 


X 

pk 


P 


p‘ 


X 


Al—k' AT—Jb—I 

= OnP + ttN-lP + . . . Ofe+ip 


( 9 ) 


whence, by subtracting both equations, we obtain the digit a^. 


The digit function, for p G N, fc G Z and x G M is then defined as [15 


dp{k,x) = 


P 


-P 


X 


k-\-l 


p 


( 10 ) 


and gives the /c-th digit of the real number x (when it is non-negative) in a positional numeral 
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system in radix p > 1. If p = 1 the digit function satisfies di(fc,a;) = di(0,x) = 0 and it does 
not relate to a positional numeral system. 


With the digit function we can express Eq. ([^ as 

Uogp hij 

X = sign(a;) E |a;|) (11) 

k=—oo 


Note that the digit function is also dehned for x negative since in that case we have, from the 
dehnition 


dp{k, -x) = p 


X 


X 



pk 


= dp(/c,p^+^ - x) 


( 12 ) 


where we have used that \_—y\ = — \y] (for any real number y). 


Most of the properties of the digit function are made apparent by plotting dp{k,x) for a hxed 
value of p as a function of x and different values for k, as shown in Fig. where ds{k,x) is 
shown for k = —2, —3, —4 in the interval x G [0,1]. The digit function is a staircase of p levels 
taking discrete integer values between 0 and p — 1. Each time that x is divisible by p^+^ the 
ascent of the staircase is broken and the level is set again to zero and a new staircase begins. 
The next proposition formalizes some of these properties. 


Proposition 1. The digit function satisfies (for n and m nonnegative integers) 


dp(fc, X + np^"*"^) = dp{k, x) 

(13) 

dp{k,p’'x) = dp{0,x) 

(14) 

dp{0,dp{k,x)) = dp{k,x) 

(15) 

,(0, n + dp(0, m)) = dp(0, n + m) 

(16) 

dp(0, ndp(0, m)) = dp(0, nm) 

(17) 

dp(n, p”") = dp(m, p*") = B {m - 

(18) 

m>n 

dp{k,n)= Y, dp(/c, j)dp(j,p'^) 

j=0 

(19) 


Proof Eq. (13) follows directly from the dehnition. We have 
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dp{k,x + = 


X + 

pk 

X 

— +np- 

pK 


-P 


X + np^^^ 


p, 


,A:+1 


X 


pK 


+ np 


-p 


X 




P‘ 


+ n 


p 


p 


X 

/c+1 


np = 


p 


p 


p‘ 


X 

:k+l 


= dp{k,x) 


( 20 ) 


Eq. (14) is also trivial to prove 


dp{k,p'"x) = 


pkx 

pk 


-p 


h 

p X 


p 


fc+1 


= \_x\ — p 


= dp(0,a;) 


Eq. (15) follows by noting that dp{k,x) is an integer 0 < dp{k,x) <p—l. Then 


dp(0,dp(fc,x)) = dp{k,x) -p 


dp{k,x) 


p 


= dp{k,x) 


( 21 ) 


since \_dp{k,x)/p\ = 0. Eq. (16) is proved by using Euclidean division, since we can always 
write m = ap + b with a, b integers and b = dp(0, m). Therefore 


dp(0, n + dp(0, m)) = dp(0, n + ap + dp(0, m)) = dp(0, n + m) 


( 22 ) 


where Eq. (13) has also been used. Eq. ([T^ is also proved in a similar way, using the distributive 


property of ordinary addition and multiplication as well. Eq. (18) is simply proved from the 


definition and Eq. (19) is a simple consequence of it. □ 


The zeroth digit of the non-negative integer n in any radix p, dp(0, n), governs the divisibility of 
n by the radix and yields the remainder of that division. In Fig. ([^ dp(0, 60) (top) and dp(0, 59) 
(bottom) are plotted vs. p. There are several interesting features that we notice immediately. 
In the case n = 60 the function touches the abscissa at the radix values that divide 60. Such 
intersections are absent in the case n = 59 (save, of course, p = 59) because 59 is prime. Of 
course, Y{p =2 dp(0. A) ^ 0 only if A is prime. 


We also see linearly decreasing branches with different slopes after a certain threshold value of 
p. This behavior is readily understood by the definition of the digit function. Note that, from 


Eq. (10), when p is sufficiently large the digit function reduces to 
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dp(0,n)=n—p if n/2 < p < n 

= n — 2p if n/3 < p < n/2 

= n —3p if n/4 < p < n/3 


etc. The function values are all contained within a triangle formed by the lines /i (p) = p — 1 
and /2 (p) = n — p. The following property follows easily 


n 


n 


n dp(o,p-1) = n(p-1) 

p=2 


{n — 1)! 


(23) 


3 Main results 


More signihcant digits exhibit even a more complex and non-monotonic behavior. Remarkably, 
though, those signihcant digits can be used to split the composition of two functions over hnite 
sets. This is the result contained in the following theorem. 


Theorem 1. (Composition-decomposition theorem.) Let g : —)■ A and f : A ^ S 

be two functions on hnite sets S consisting of all integers G [0,p — 1] and A consisting of all 
integers G [0,p^ — 1] with p, iV G N and let k G . We have 


f{g{k)) = dp 


\ 

9{k), Y. 

n=0 J 




p; 




p 


EU P"/(") 

p9(^)+l 


(24) 


Proof. From Eq. (14) we have 


dp 


\ 

9{k), Y 

n=0 J 


0 , 


p ^-1 

pn-g{k) 

n=0 


f{n) 


(25) 


This function extracts the zeroth digit of the quantity I]n=o^p”'“^^^^/(n), i.e. the value of /(n) 
for which the exponent of the accompanying power of p within the sum is zero. But this happens 
when n = g{k) from which the result follows. □ 


We have now all what we need to prove our main result. 









Theorem 2. (Universal map for CA in closed form.) Eq. 0 is equivalent to 


X 


t+i 


= dp ^ = 


L k=—r 


R 




-P 


R 


i+u! p'^+^xR'" 


( 26 ) 


Proof. Let us take N = I + r + 1. Eq. Q can be seen as the composition of two functions. The 
Erst one g{k) = sends the N integers xi'^^ (each G [0,p — 1]) to one integer 

nl G [0,p^ ~ 1] = — 1] and hence qualihes as a valid function g{k) for the composition- 

decomposition theorem. The second one, f{n) = a„ sends the integers n G [0,p^ — 1] to one 
integer G [0,p — 1]. Thus, by the composition decomposition theorem we have, by replacing 
g{k) and f{n) under this identification 


X 


t+i 


l pl+r+l 

fc+r i+fc „n. 


= dp , 51 P"‘ 


ik=—r 


n=0 


from which, by using Eq. (|^, the result follows. □ 


(27) 


The above result replaces a sum over — 1 terms by just only one term evaluated on 

the neighborhood of the site at location j. Furthermore, we do not need to give any table of 
configurations as input but just only the Wolfram code (a nonnegative integer R G 
and the neighborhood parameters p, I and r directly. For example, the celebrated Wolfram’s rule 
110 takes the simple form 


1 


X 


j _ 


= I 5: 2^+^a;r5ll0 

k=-l 


(28) 


The maps for all Wolfram’s 256 CA with / = r = 1 and p = 2 are thus simply obtained as 


1 


X 


J _ 


*+, = d2| 5: 2^+^a;r,i? 

k=-l 


(29) 


by giving an input R G [0, 255]. 


As a further illustration, the rule = d 2 1771466136^ is plotted in Fig. 

1^ starting from an arbitrary initial condition. This rule displays many coherent soliton-like 
structures that are created and annihilated against a quiescent background. 


With Eq. (26) alone we are able to explore the whole CA rule space within computational 
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reach and go even beyond, because the mathematical expression is valid for arbitrarily large 
parameter values regardless of our computational power to carry out the simulations of a given 
rule. 
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1 The i3-function for real-valued arguments (top) and for ?/ = ±e constant in the 

second argument (bottom). [13] 

2 The digit function d^^kyx) plotted for k = —2, —3, —4 (from top to bottom) in 

the interval x G [0,1] HH 
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4 Spatiotemporal evolution of the rule x{^i = d 2 ^ 1771466136^ 

starting from an arbitrary initial condition on a ring of 400 sites and 200 time 
steps. Time flows from top to bottom. |T6] 
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Fig. 1. The i3-function for real-valued arguments (top) and for y = ±e constant in the second argument 
(bottom). 
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Fig. 2. The digit function d 3 {k,x) plotted for k = —2, —3, —4 (from top to bottom) in the interval 

X G [0,1] 
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Fig. 3. The digit function dp(0,60) (top) and dp(0,59) (bottom) plotted vs. p. 
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Fig. 4. Spatiotemporal evolution of the rule = d 2 1771466136^ starting from 

an arbitrary initial condition on a ring of 400 sites and 200 time steps. Time flows from top to bottom. 
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